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Molecular clouds, the birthplaces of stars in galaxies, form dynamically from the diffuse atomic 
gas of the interstellar medium (ISM). The ISM is also threaded by magnetic fields which have a 
large impact on its dynamics. In particular, star forming regions must be magnetically supercrit¬ 
ical in order to accomodate gas clumps which can collapse under their own weight. Based on 
a parameter study of three dimensional magneto-hydrodyamical (MHD) simulations, we show 
that the long-standing problem of how such supercritical regions are generated is still an open 
issue. 


1 Introduction 

Present day stars form within the densest regions of molecular clouds (MCs) and giant 
molecular clouds (GMCs), in gravitationally unstable cores and clumps. Our common un¬ 
derstanding is that those MCs and GMCs form from the diffuse, atomic (HI) gas within 
timescales of less than lOMyiEEl The generation of filaments and substructures within 
GMCs is primarily controlled by magnetic fields and turb ulencd^l^ In particular, mag¬ 
netic fields are an elemental part of the interstellar mediurt El ^l^ l ^l*°P^ ] which have a large 
impact on the dynamics of the ISM on various spatial scaleJi^IIIIlIIIl as the magnetic en¬ 
ergy density is comparable to the thermal energy density of the ISIVP. 

One long-standing issue is the formation of supercritical clumps and cores. Sim¬ 
ilarly to thermal pressure, magnetic fields prevent contraction of otherwise (thermally) 
self-gravitating gas clumps if the magnetic fields are strong enough. Therefore, gaseous 
overdensities must be magnetically supercritical, quantified by the mass-to-flux ratio, p, 
to collapse and to subsequently allow the formation of stars. 

Already in 1956, Mestel & Spitzer realised that molecular clouds should be magneti¬ 
cally subcritical assuming field strengths that correspond to the equipartition of magnetic 
and kinetic energy density within the ISIVP^I. To generate supercritical cloud cores out of 
those subcritical conditions, they suggested that the non-perfect coupling between charged 
particles and neutrals, i.e. the ambipolar diffusion (AD) drift, could locally reduce the 
mass-to-flux ratio which allows the cloud to break up and to form stars . Fo r a long time 
this was the standard theory of star formation out of the magnetised isiviim In this fairly 
static “standard model” of magnetically-supported, AD-mediated supercritical cores, low- 
mass stars would form by the the slow gravitational contraction of isolated cores containing 
a very small fraction of the clouds’ mass. This picture would also account for th e very low 
observed global star formation efficiency (SFE) of giant molecular cloud^SEHl The slow 
contraction results from the typical timescale for ambipolar diffusion, ^ad. which is an 
order of magnitude larger than the free-fall time, tg, of individual cloud cores (their ra¬ 
tio is about fAD/iff ~ 10 (xe/10”^), where is the ionisation fraction). On the other 
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Figure 1. These observational data summarise the main motivation for our proposed study: How do subcritical 
(HI) clouds become supercritical (H 2 ) clouds? Our previous studies have shown that it is everything but trivial 
to build up supercritical clouds out of the magnetised interstellar medium, because the mass-to-flux ratio is fairly 
well conserved, even in the presence of ambipolar diffusion and enhanced non-ideal MHD setupJSES [From 
Crutcher (20I2pl| 


hand, rather recent observations by Crutcher (2009^ of individual cloud cores includ¬ 
ing Zeeman measurements to determine their magnetic field distribution indicate that ide¬ 
alised models of ambipolar-diffusion driven star formation are unlikely to be operative. In 
idealised models with ordered background magnetic fields, efficient ambipolar diffusion 
would lead to a local increase of the mass-to-flux ratio towards the centre of cloud cores 
which is not seen in their observed sampldSl (but see also Bertram et al. (2012jSlon the 
difficulty to interpret those observations). 

However, present-day models of star formation also account for th e fact th at molecular 
clouds are also pervaded by supersonic random motions, i.e. turbulencd^lEEIl Eventually, 
this resulted in a paradigm shift of the theory of star formation where magnetic fields only 
play a minor role and supersonic, super-Alfv enic turbulence controls the star formation 
efficiency within molecular cloudJlElIlZlllIlH As a consequence, the magnetic fields are 
expected to be highly disordered rather than being an ordere d bac kground field. Hence, 
idealised models of ambipolar diffusion drift should not appl>l^^l^. Additionally, the AD 
characteristic timescale is expected to decrease in t his ca sj ^^ l ^^l^"^l^^ land other diffusive ef¬ 
fects like turbulent reconnection might be o perativj^^l^ In deed, a number of studies have 
suggested that both MC and their clumpglHE^lIElElEll are close to being magnetically 
critical, with a moderate preference for being supercritical. Moreover, recent compila¬ 
tions of observational data show that cloud cores and clumps with column densities of 
N ^2 X 10^^ cm“^ are essentially all supercritical (see Fig.[^. 

Whether those supercritical cloud cores and clumps are the result of ambipolar dif¬ 
fusion together with random motions in the ISM is far from being certain and has to be 
investigated further. For instance, recently Heitsch & Hartmann (2014j^ argued in their 
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Figure 2. Setup of the initial conditions. Left: The well studied head-on collision where both flows ai‘e in the 
direction of the background magnetic field. Right. Flows with an inclination angle with respect to the background 
magnetic field. This setup resembles the impact of external driving of such flows (e.g. driving by SN blast waves). 
[From Kortgen & Banerjee (2015pll| 


parameter study that ambipolar diffusion in concert with turbulence is unlikely to con¬ 
trol the formation of supercritical cores and hence star formation. They again propose an 
alternative scenario where large scale flows are the main driver to generate supercritical 
cores. This idea, where supercritical clouds could be assembled from large scale flows was 
already discussed in Mestel & Spitzer (1956)21 as an alternative to the AD-mediated sce¬ 
nario and to avoid the “magnetic flux problem”. But only in combination with supersonic 
turbulence this scenario becomes more feasible because gravitational frag mentat ion could 
be suppressed during the assembly of the clouds by those turbulent motionJ^IEEl This ac¬ 
cumulation idea would also support a number of recent observations which show that mag- 
ne tic fields are dy namically important on all scales in the Milky Way and other spiral galax- 
iej 50 | i 2 |i 4 | 5 i | i 3 |i 5 [ jg particular evident from Fig. 1: The low column density HI gas 
is magnetically subcni\c?\, whereas clouds which exceed columns of iV > 2 x 10^^ cm“^ 
are magnetically supercritical. 

In the presented numerical parameter study, we investigated the possibility of diffusion 
mediated generation of supercritical clouds showing that it is unlikely that such unstable 
clouds can be build up from subcritical Hl-clouds. 


2 Numerical method and initial conditions 

For these studies we used the FLASH adaptive mesh refinement (AMR) codj^^l ^ In ad- 
dition to the basic ideal MHD equations (for which we employ the Bouchut solvetEHHHj) 
we also used the ambipolar diffusion module developed from Duffin & Pudritz (2008)21 
Additionally, self-gravity as well as heating and cooling processes were included in those 
simulations. For the latter, we followed the treatment by Koya ma & Inutsuka (2002)21 
(an analytic simplification of their detailed calculation ii j59 | 60 | 6i j) capture the build-up 
of self-gravitating cores within the molecular clouds we used sink particle^ in addition 
to the Jeans refinement criterion (i.e. the Truelove criterior|21). In particular the detailed 
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sink particle approach allows us to unambiguously identify supercritical, collapsing re¬ 
gions which are importent for our studies quantifying the star formation ability from the 
magnetised ISM. _ 

Our initial setups for those studies are similar to the ones described irl^MIEl] (see also 
Fig. 1^ where the build-up of molecular clouds is modeled by the collision of cylindrical 
streams of warm neutral HI gas (WNM). Each flow is ^ = 112 pc long and has a radius of 
r = 64 pc. The bulk flows are slightly supersonic with typical Mach numbers of Alf = 2. 
On top of those bulk motions, a turbulent velocity field is superimposed which triggers 
initial instabilities like the non-linear thin-shell instability (NTSlj^and subsequently leads 
to fragmentation of the cloud. The initially uniform magnetic field has a strength of B = 
{3,4, 5} /rG corresponding to mass-to-flux ratios of /t//icrit ~ 1,0.7, 0.6 if the critical 
value /Tcrit ~ 0.13/v^is applied^. 

Furthermore, we also studied the impact of an oblique angle of the flows with respect 
to the background magnetic field (see right panel of Fig. |^. Those oblique flows are more 
realistic than the head-on flows and could be generated, for instance, by supernova shock 
waves and by the gravitational potential of spiral arms. The motion of the flow at an in¬ 
clination with resp ect to th e magnetic held results in enhanced magnetic diffusivity (by 
numerical diffusiorI^3^3Ilj). Again, those flows resemble streams of the WNM in a ther¬ 
mally bistable conflguratioiP^. The flows are studied with different oblique angles, which 
are varied from 10° to 60°, different initial magnetic fields strengths and different strength 
of the initial turbulence ranging from subsonic to supersonic velocity fluctuations. For 
details on the numerical setup and our initial conditions see Kortgen & Banerjee (2015j^. 

3 Results 

As can be seen from Fig.[^ the different initial held strengths have significant implications 
for the resulting dynamical behaviour of the molecular cloud. The main difference comes 
about in efficiency to form stars (or not). In the case of a rather weak background held of 
3 fjG supercritical star forming clumps can be generated whereas in the case of a slightly 
stronger, but more realistic, magnetic held star formation is fully suppressed. Note that, 
due to the oblique flows with an angle of 60° the effective mass-to-flux ratios are 0.73 
in the 3 pG case and 0.44 in the 5 /rG case. That means that both cases are initially sub- 
critical, but only in the cases of the weak magnetic held locally supercritical clumps are 
assembled due to sufficient flux loss. 

An interesting point is also the field morphology. In the weak magnetic field case 
the field structure in the dense regions is clearly separated from the large scale magnetic 
field, whereas in the strong field case the field morphology is almost unaffected compared 
to the initial configuration (see the blue stream lines of Fig. [^. From an observational 
point of view, the field structure and its dynamical importance within molecular clouds 
is still debated. On the one hand, some multi-scale polarisation data indicate that mag¬ 
netic fiel ds in GMC s are essentially just dragged in from larger scales and are dynamically 
importanCnnUmiD On the other hand, Zeeman measurements of individual cloud cores 
together with analyses of nume rical simulations indicate rather weak fields that might not 
be dynamically import an i^^ l ^H With our subsequent studies on cloud formation on cloud 
scales including a more detailed modelling of ambipolar diffusion we hope to clarify this 
issue. 
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Figure 3. Results from colliding flow simulations investigating the formation of molecular clouds. Here, the 
flows collide with an oblique angle of 60°. Left panel: The weak field case (3 pG). The Right panel shows the 
same situations in the case of a stronger background magnetic field of 5 pG. In the case of a weak magnetic 
field supercritical cloud cores can form that allow the formation of stars (marked with black dots). Stronger 
initial magnetic fields prohibit the formation of stars even in the case of large oblique angles of the flows. The 
blue stream-lines indicate the magnetic field morphology in the projected 2D plane. [From Kortgen & Banerjee 
(2015pil| 

In Fig. 1^ we quantify the main results by means of histograms in the N-B-p\ane from 
our colliding flow studies for various initial conditions. Only in the initially marginally 
subcritical case (B = 3 fiG) we observe signs of star formation within supercritical core]^ 
For slightly stronger initial magnetic fields {B^A /iG) no supercritical cloud cores are gen¬ 
erated, hence there is no star formation activity, regardless whether ambipolar diffusion is 
active or the flows collide with an oblique angle. Nevertheless the results of those simula¬ 
tion show the observed behaviour in the low column regime (N < 10^^ cm“^), where gas 
assembles along field lines without changing the field strength by much (see also the latest 
analysis from PLANCK observations of individual molecular cloudsOnly within su¬ 
percritical, self-gravitating cores the magnetic field gets enhanced by compression due to 
flux freezing. Furthermore, we observe that star formation is immediately initiated, when 
the gas becomes supercritical promoting a picture of “rapid” star formatioJ^. 

4 Conclusions 

Here, we summarised our recent results from MHD simulations of colliding flows with 
varying initial conditions on the possible formation of supercritical cloud cores from sub- 
critical initial conditions. Although, dense clouds are easily formed within colliding flow 
scenarios due to thermal instability, the generation of supercritical clumps are largely de¬ 
termined by the initial conditions. Furthermore, increasing initial turbulence lead to lower 
masses of the cores and clumps because the HI streams become less coherent. Otherwise, 
increasing magnetic field strength lead to more massive molecular clouds, which never¬ 
theless do not become supercritical. Oblique flows still lead to cloud cores with masses 

®If we assume ~ 0.13/\/G for spherical core m we get /J.//.icrit = 0.97 (3 pG/B) for our head-on 
colliding flow configurations. 
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Figure 4. Results from colliding flow simulations with various different initial conditions. Shown are histograms 
of the line-of-sight field strength -Blqs function of the column density N. From left to right: B = 3 fiG,B = 
4/iG and B = b fiG, respectively. Top: $ = 0°, Bo ttom : $ = 60°. Different line colours denote different 
times. Also shown ai'e the criticality condition (red lind^^FI), corrected for projection effects (black linJ^, and 
assuming equipartition of turbulent and magnetic fields (blue lin£^. Colour coded is the mass distribution within 
this two parameter space, (from Kortgen & Banerjee (2015pl| 


comparable to what has been observed recently. But starting with subcritical HI flows, 
in no case the magnetic flux loss is sufficient to allow the build-up of supercritical cloud 
cores. Generally, increasing inclination of the flows lead to increasing diffusivity of the 
magnetic field. Again, regardless of the variation of the inclination, no tendency for faster 
accumulation of gas or faster transition to thermally dominated regions was seen in our 
simulations. 

We therefore stress the role of magnetic fields in the context of molecular cloud and 
star formation. We point out the complete lack of supercritical regions for realistic initial 
field strengths. From the observational side, HI clouds may be supercritical as a whole, but 
their observed, dense subregions be subcritical. 

Hence, the question remains, how magnetically supercritcal cloud cores are formed? 
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